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(57) Abstract 

A method and apparatus are disclosed for determining more accurately, in a radio environment with time dispersion, the distance 
between a radio receiver (100) and radio transmitter (130), by special processing of received radio signals that have been transmitted 
repeatedly from the same radio transmitter (130) and are possibly subject to muKipath propagation. The Time of Arrival (TOA) of the 
received radio signals is repeatedly estimated using channel power profiles (step 402). A TOA value near the minimum occurring TOA is 
selected (step 402) wherein each estimated TOA is derived from incoherent integration of a randomly chosen number of the received bursts 
having the same known bit sequence (step 404), in order to eliminate the influence of noise. 
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METHOD AND SYSTEM FOR DETERMINING 



POSITION OF MOBILE RADIO TERMINALS 



-> BACKGROUND OF THE INVENTION 
Technical Field of the Invention 

The present invention relates in general to the mobile radio communications 
field and, in particular, to a method and system for use in determining the 
geographical position of a mobile radio terminal. 

, Description of Related Art 

Commonly-assigned U.S. Patent Application Serial No. 08/894,466 to 
Hagerman et al., describes a method for determining the position of a mobile radio 
terminal (hereinafter mobile station or MS), which makes use of uplink Time of 
Arrival (TOA) and Direction of Arrival (DOA) measurements. Slight mention is 
made in that disclosure about repeating the same uplink signal a number of times, 
in order to improve the accuracy of the TOA determinations. However, the 
disclosure does not provide any details about how this improvement can be 
achieved. 

Commonly-assigned U.S. Patent Application Serial No. (Attorney docket 
no. 27946-00353), which has been filed on the same day as the present application, 
discloses a method and apparatus for providing a "prime" uplink signal from a MS 
whose position is to be determined. The present invention is related by subject 
matter to the above-described Patent Applications and teaches how to process a 
repeated uplink signal, such as, for example, the "prime" uplink signal from a MS 
whose position is to be determined. 

The CME 20 System, which is a Global System for Mobile 
Communications (GSM) implementation manufactured by Telefonaktiebolaget LM 
Ericsson, sends uplink messages with a Timing Advance (TA) to compensate for 
the radio signal propagation delay caused by the distance between the radio base 
station (BS) and the MS. The TA value is determined by a procedure of selecting 
the smallest TA value from a number of measurements made by the BS on several 
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access bursts that have been transmitted by the MS. For TA value determinations, 
this procedure solves the radio signal multipath propagation problem (i.e., when the 
signal propagates along a line of sight path and also along one or more reflected 
paths). This procedure is hereinafter referred to as the "minimum time" method. 
5 SUMMARY OF THE INVENTION 

' A problem addressed by the present invention is how to process a number 
. of identical digital signals transmitted from a radio transmitter (TX), which are 

subject to noise and multipath propagation when received in a radio receiver (RX), 
in order to more accurately determine the distance between the TX and RX than 
1 o what is possible using the above-described "minimum time" method. 

Another problem addressed by the present invention is how to process the 
above-described digital signals if the signal parameter being determined in the RX 

is the TOA. ■ x 

Still another problem addressed by the present invention is how to process 
15 the above-described digital signals if the signal parameter being determined in the 
RX is the TOA, and there are TOA and DOA measurements being made in a BS. 

It is, therefore, an object of the present invention to provide a method and 
apparatus to more accurately determine the distance. between an RX and TX, by 
processing radio signals that have been subject to noise and multipath propagation, 
2 o wherein the radio signals are being transmitted repeatedly from the same TX and 

received in an RX. 

It is another object of the present invention to provide a method and 
apparatus in accordance with the above-described object, whereby the signal 
parameter being processed is TOA. 

2 5 it is yet another object of the present invention to provide a method and 

apparatus in accordance with the above-described objects, whereby the signal 
parameters being processed are TOA, TDOA and DOA. 

The foregoing and other objects are achieved by a method and apparatus for 
improving time measurement accuracy, which includes transmitting "M" samples 

3 o of the same signal from a TX, receiving the "M" signal samples along with 

multipath components and noise, for each of the received "M" samples, 
determining in the RX an estimated channel power profile (CPPi), selecting "N" 
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out of the "M" received samples, performing incoherent integration of the CPPi for 
the "N" samples, which results in an integrated signal, ICPP(Ni), determining if the 
signal-to-noise quality of the ICPP(Ni) is greater than or equal to a predetermined 
threshold value, and if not, improving the signal-to-noise quality of ICPP(Ni) as 
required, by redoing the incoherent integration with successively one additional 
received sample CPPi until the signal-to-noise quality of the ICPP(Ni) is greater 
than or equal to the predetermined threshold value, determining the TOA(i), 
including the case of determining TOA(i) from the maximum signal amplitude, and 
entering that TOA(i) value into a diagram that shows a frequency of occurrence as 
a function of TOA(i), repeating the whole procedure "X" times by selecting a new 
combination of "N" out of "M" samples, which results in "X" additional points in 
the frequency of occurrence diagram, and reading the minimum value TOA(min) as 
the time value having "z" of all occurrences with higher TOA(i) values and "1-z" 
of all occurrences with lower TOA(i) values, where z >0.7. 

An important technical advantage of the present invention is that the 
method for improving the quality of the signals, which requires the addition of 
many signals (incoherent integration), is performed in such a way that also allows 
the minimum TOA to be determined as in the "minimum time" approach described 
above. 

Another important technical advantage of the present invention is that the 
elimination of noise thus provided allows more remote BSs to take part in 
determining the position of a MS, which not only improves the accuracy of each 
measurement made, but also improves the accuracy overall because of the 
increased number of BSs that can participate in the position determination process. 
BRIEF DESCRIPTION OF THE DRAWINGS 

A more complete understanding of the method and apparatus of the present 
invention may be had by reference to the following detailed description when taken 
in conjunction with the accompanying drawings wherein: 

FIGURE 1 is a schematic block diagram of a radio BS that can be used to 
perform measurements for determining TOA, in accordance with a preferred 
embodiment of the present invention; 

FIGURES 2A-2E illustrate a series of wave-forms associated with a method 
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for determining a sample TOA(i), which can be implemented in accordance with 
the preferred embodiment of the present invention; 

FIGURE 3 is a diagram that can be used to derive an estimated minimum 
Time of Arrival or TOAe based on seventy bursts from a TX, in accordance with an 
5 exemplary embodiment of the present invention; 

FIGURE 4 is a flow diagram that can be used to implement an exemplary 
embodiment of the present invention; 

FIGURE 5 is a diagram of an alternative representation for the diagram of 
FIGURE 3, and which shows a cumulative distribution curve for the same bursts; 
10 FIGURE 6 is a flow diagram that illustrates a first auxiliary method, for 

rejecting interference in accordance with the present invention; and 

FIGURE 7 is a flow diagram that illustrates a second auxiliary method, for 
detecting the measurement region in the time domain, in accordance with the 
present invention. 
1 5 DETAILED DESCRIPTION OF THE. DRAWINGS 

The preferred embodiment of the present invention and its advantages are best 
understood by referring to FIGURES 1-7 of the drawings, like numerals being used for 
like and corresponding parts of the various drawings. 

Essentially, the present invention provides a method for making time 
2 0 measurements which can be used in determining the geographical position of a mobile 
radio terminal. As such, the time measurements involved can be TOA measurements 
or Time Difference of Arrival (TDOA) measurements. The following description 
discloses an embodiment of the invention wherein TOA measurements are preferably 
used. 

2 5 Specifically, FIGURE 1 is a schematic block diagram of a radio BS that can 

be used to perform measurements for determining the TOA of signals, in accordance 
with a preferred embodiment of the present invention. As shown, the BS 100 can be 
connected for communication with a wire-bound network (not explicitly shown) via 
a mobile services switching center (MSC) 125. As shown, the BS 1 00 performs TOA 

3 o measurements for a mobile radio terminal (MS) 1 30 over a radio connection, of which 

the two uplink paths 101, 102 are shown that provide diversity. For a cellular mobile 
radio BS, it is common to include two receiving antennas (e.g., 103 and 104) which 
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provide this uplink diversity. Associated with each antenna 103 and 104 is a 
respective receiver 111a and 111b, radio frequency (RF) demodulator 1 12a and 1 12b, 
and intermediate frequency (IF) demodulator 113a and 1 13b. A received frequency 
synthesizer 114 is connected to each of the RF demodulators 112a and 112b and 
5 control unit 115. The received frequency synthesizer 114 enables reception of 
particular radio channels using these connections. Demodulated signals from the IF 
demodulators 113a and 113b are coupled to an equalizer 135 over respective 
connections 116 and 117. Although not explicitly shown, the equalizer is also 
connected to other components in the BS 100. 

10 The demodulated signals from the IF demodulators 1 13a and 1 13b are also 

coupled to a modified-RX (ModRX) 118. The ModRX 118, the function and 
operation of which are described in the above-mentioned U.S. Patent Application 
Serial No. 08/894,466 to Hagerman et al., functions to measure the TOA of the signals 
received from the MS 130, and reports the results of these measurements to the MSC 

15 125 via a connection to the control unit 115. Preferably, the processing of the TOA 
measurements is performed by the processor 1 19 in the ModRX 118. Alternatively,, 
the processing of the TOA measurements can be performed by a processor in control 
unit 1 15, or in the wire-bound network (e.g., in the MSC 125). The equalizer 135 
combines the two demodulated signals on connections 116 and 1 17, in order to take 

2 0 advantage of both the incoming line-of-sight and reflected radio energy. On the other 
hand, the ModRX 118 processes these two demodulated signals separately, because 
signal fading is a desired effect when measuring TOA. In other words, the fading 
eliminates reflected beams for some received bursts, and thus provides only the 
desired line-of-sight signals. 

2 5 FIGURES 2A-2E illustrate a series of wave- forms associated with a method for 

determining a sample TOA(i), which can be implemented in accordance with the 
preferred embodiment of the present invention. The wave-forms shown in FIGURE 
2 A represent one burst of the received signal on connection 116 (or 117), which 
contains variable data 1 and data 2, as well as a training sequence (known in advance) 

3 0 in the middle of the burst. The solid wave- form represents a direct (line-of-sight) 

received signal, and the dotted wave- form represents a received reflection of the burst 
(thus slightly delayed in time). 
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The wave-forms shown in FIGURE 2B illustrate the theoretical responses of 
a correlation between the incoming direct and reflected bursts of the signal on 
connection 1 16 (or 1 17), and the known training sequence available in the ModRX 
1 1 8 for determining TOA. The combined correlation response ("Channel Impulse 
Response" or CIR) of these bursts as processed in the ModRX 118, is shown in 
FIGURE 2C. The suffix "i" indicates a particular received burst. 

The CIRi wave-form shown in FIGURE 2C is a complex value comprising 
amplitude and phase. The central part of the CIRi is given by a limited number of 
samples (e.g., 5), but a complete signal as shown in FIGURE 2C is reconstructed by 
the processor 1 19 using a conventional interpolation technique. The wave-form 
shown in FIGURE 2D illustrates the squared absolute value of the CIRi, which is 
referred to as the "estimated channel power profile" or CPPi. The wave-form depicted 
for the CPPi is optimistic, and in actual cases often contains a large noise component. 
Consequently, N out of M (e.g., M=70) CPPi are integrated into a signal ICPP(Ni), as 
shown in FIGURE 2E. 

The integrated signal ICPP(Ni) has a higher peak value, h, for the desired 
signal, and a bias value, b, for the undesired noise, but with a spread sigma for the 
noise that is not amplified to the same extent by the integration. The bias value, b, is 
eliminated during processing, and the quality "Q" of the ICPP(Ni) is defined as 
h/sigma. A suffix, Ni, assigned to the Integrated Channel Power Profile ICPP(Ni) 
indicates that a certain selection "i" of N bursts has been integrated. A TOA(i) value 
is associated with the ICPP(Ni) and is entered as one of 100 points, for example, in the 
diagram shown in FIGURE 3. 

FIGURE 3 is a diagram that can be used to derive the estimated minimum time 
of arrival (TOAe) based on 70 bursts, in accordance with an exemplary embodiment 
of the present invention. Although FIGURE 3 discloses a diagram that can be used 
to plot TOA(i) versus frequency of occurrence, it is shown as a diagram for illustrative 
purposes only and is not intended to so limit the present invention. In an automated 
processing environment, the diagram of FIGURE 3 can be implemented, for example, 
as a multi-dimensional matrix or array, with the rows (or columns) including each 
TOA(i), and the columns (or rows) including each occurrence. An alternative 
representation is given in FIGURE 5, which shows a cumulative distribution curve for 
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the same 70 bursts. 

The present invention preferably uses M (e.g., 70) received bursts containing 
at least a part that has identical known contents of bits in all received bursts (e.g., a 
training sequence). However, note that a different embodiment can use (e.g., 70) 
received bursts containing different (non-identical) information, but information still 
known to the receiver, which enables the processor to produce CIRi by correlating 
each burst with its known contents. 

Using a first procedure, the noise is reduced to a preselected level by "adding" 
(e.g., incoherent integration of the CPPi) a required number of N bursts which are 
randomly selected from among the M received bursts, and determining the TOA(i) for 
this sum, ICPP(Ni), of N bursts. This first procedure is repeated a plurality (e.g., 100) 
times, which for this example, results in 100 values of TOA(i) which are entered into 
the diagram show in FIGURE 3. The primary advantage of this first procedure is that 
the random selections, Ni, will sometimes comprise bursts that contain mostly (or 
only) bursts received on the shortest path (line-of-sight path). Using a second 
procedure, the minimum TOAe is read in the diagram shown in FIGURE 3. However, 
the shortest values of TOA(i) are excluded because they may be due to noise and 
multipath effects displacing peaks. Empirical information shows that a suitable 
reading for TOAe is where 5% of the occurring TOA(i) are earlier than the TOAe, and 
95% of the occurring TOA(i) are later than the TOAe. 

FIGURE 4 is a flow diagram that illustrates an exemplary method 400 that can 
be used to implement the present invention. At step 401, M bursts of the received 
signal at connection 1 16 (or 117) that contain a known bit sequence are measured by 
ModRX 118. For example, these bursts can be access bursts sent on the uplink from 
the MS 130 as a response to an intra-cell handover order, or a handover order back to 
the pertinent (old) traffic channel, with the order coming from a serving BS (e.g., 100) 
in a GSM network. If these access bursts are not answered by the BS 100, up to 70 
bursts will be transmitted by the MS 130 before a time-out occurs and places the MS 
back on the previously used (old) traffic channel. The demodulated signals on 
connections 116 and 117 both include the received access bursts, or together 2 x 70 
= 140 bursts maximum, each of which can be treated as an independent burst with 
respect to the method of the present invention. 
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At step 402, the received channel impulse responses, CIRi (as shown in 
FIGURE 2C), which are each associated with one particular burst and contain complex 
information, are each processed to form an estimated channel power profile, CPPi (as 
shown in FIGURE 2D) equal to the squared absolute value of the CIRi. Consequently, 
there will be M values of CPPi. At step 403, the processor 1 1 9 randomly selects N out 
of the M values (wave-forms) of the CPPi starting with N=10, for example. 

At step 404, an incoherent integration of the N=10 selected wave- forms CPPi 
results in one integrated channel power profile ICPP(Ni), as shown in FIGURE 2E. 
At step 405, a check is made to determine whether the noise has been reduced enough 
for the derived ICPP(Ni) to have a minimum quality defined by the ratio h/sigma 
(compare FIGURE 2E). If the h/sigma ratio is less than a predetermined threshold, the 
derived ICPP(Ni) has to be improved by repeating the steps 403-405 using one 
additional CPPi randomly selected from among the not yet selected CPPi. At step 
407, the minimum quality has been reached, and the resulting TOA(i) (as shown in 
FIGURE 2E) is entered into the diagram shown in FIGURE 3. The flow then 
continues to step 408. 

At step 408, the question is answered about whether or not the predetermined 
number of TOAi values has been reached. If not, the steps 403-407 are repeated as 
described above, until the required predetermined number of TOAi values has been 
reached. 

If so, at step 409, the lowest TOA(i) in the diagram of FIGURE 3 is selected, 
because this value of TOA represents the TOA of the burst as received on the shortest 
path (line-of-sight path). However, the selection is performed according to the 
description for FIGURE 3 or 5 above, thereby rejecting the influence of noise. 

In another embodiment of the present invention, TDOA is measured instead 
of TOA. For the TDOA measurements, the procedure is the same as the procedure 
described above for TOA, except that time measurements are made relative to another 
signal rather than relative to a time reference. If TDOA measurements are performed, 
the cross-correlation of two received signals can be performed. As such, the TDOA 
estimate can be the highest peak of the cross-correlation output. 

Also, if TDOA measurements are performed, the above-described inventive 
incoherent integration approach can also be used. However, two base stations should 
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be involved, in order to have two signals to perform the TDOA measurement. As 
such, the cross-correlation output can be substituted for the CIR, and the cross- 
correlation power profile can be substituted for the CPP (where the CPP represents the 
squared absolute value of the cross-correlation output). 

In another embodiment, if DOA measurements are performed, a spatial 
signature can be estimated from, for example, the received signal from a directibnal- 
or narrow-beam antenna (e.g., an adaptive antenna or fixed lobe antenna. As such, the 
DOA* estimate can be the highest peak of the estimated spatial signature. 

The above-described inventive incoherent integration approach can also be 
used if DOA measurements are performed. As such, the spatial signature can be 
substituted for the CIR, and the spatial signal power profile can be substituted for the 
CPP (where the CPP represents the squared absolute value of the spatial signature). 

In both the TDOA and DOA measurement cases, the line-of-sight (LOS) signal 
component corresponds to the "center of gravity" of the ICPP (i.e., the LOS 
component does not correspond to the lower part of the correlation output, such as in 
the TOA measurement case). Consequently, there would be little or no advantage in 
performing incoherent integration with less than the total number, M, of the received 
bursts. Therefore, in these cases, all bursts are used in the incoherent integration 
process, and the TDOA and DOA estimates are the "center of gravity" of the ICPP. 



In yet another embodiment of the present invention, an antenna array is used 



are weighted by their quality (i.e., by their h/sigma value), which is an indication of 
the noise involved. Consequently, remote BSs receiving weak signals are given less 
influence over a position determination process. 

FIGURE 6 is a flow diagram that illustrates a first auxiliary method, for 
rejecting interference in accordance with the present invention. This method teaches 
how to overcome disturbances from other mobile terminals (e.g., whose position is not 
being determined). For example, in the GSM, there are eight different training 
sequences used which have a certain correlation to each other, and can cause 



in a BS for measuring beams from different directions separately. Subsequently, the 
line-of-sight beam is selected for the TOA measurements (shortest TOAi). 

In still another embodiment of the present invention, the occurrences of TOA(i) 
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disturbances to the measurements performed for a mobile terminal whose position is 
being determined. 

At step 601 of the method, the channel impulse responses, CIRj, are estimated 
for all possible interferer training sequences on all M received bursts. At step 602, the 
5 corresponding channel power profiles, CPPj, are formed, by taking the squared 
absolute value of the CIRj . At step 603, the average carrier-to-interference ratio, C/I, 
is estimated for each such interferer training sequence, while treating the "C" part as 
the power of the interferer of interest and the "I" part as the power of other signals. 
At step 604, a determination is made about whether the disturbance is to be rejected 

0 (i.e., subtracted) from the received bursts, which is the case if the C/I exceeds a first 
threshold value, Tl . At step 605, those training sequences that meet the condition of 
step 604 are convoluted (symbolized by a "*") with the CIRj, and subtracted from the 
received bursts. 

FIGURE 7 is a flow diagram that illustrates a second auxiliary method, for 
5 detecting the measurement region in the time domain, in accordance with the present 
invention. For example, this method is particularly useful in cases of high noise 
compared to the signal of interest. As such, before searching for correlation peaks (as 
shown and described above with respect to FIGURE 2B), there can be a problem with 
determining the approximate time when the training sequences are found among the 

1 o noise. This second auxiliary method comprises incoherently integrating all M bursts 

received, which results in a relatively wide pulse that emerges from the noise, and 
indicates the approximate time when the training sequence of the measured mobile 
terminal (whose position is being determined) is to be found. 

At step 701 of the method, the channel power profiles, CPPi, of the received 

2 5 carrier are formed for all M bursts. At step 702, the incoherent integration of all the 

CPPi is performed, and the highest peak corresponds to the location of the training 
sequence. At step 703, the C/I is estimated, while treating the "C" part as the power 
of the carrier, and the 'T part as the power of other signals. At step 704, a 
determination is made about whether the result from steps 701-702 is acceptable 

3 0 (C/I>T2), or the noise is strong enough such that additional steps have to be taken for 

an acceptable result. If, at step 704, the value of C/I is above the threshold T2, the 
result from steps 701-702 is accepted. Otherwise, the bursts with high influence from 
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noise are weighted down as follows. At step 705, the bursts with high energy contents 
(which are likely due to noise) are discriminated by assigning them a low weight 
factor, Wi\ and the bursts with low energy contents are assigned a high weight factor 
for Wi\ At step 706, a sum of all M weighted bursts is formed in accordance with the 
formula: 

m . . 

H(Wt'*CPPQ. (1) 
i=i 

The highest peak of this sum corresponds to the location of the training sequence. 

The principle of assigning weights to the received bursts depending on the C/I 
can also be applied to the processing of measured bursts in order to determine a time 
of arrival, TOAe. The training sequence can be used to estimate the power of the 
signal as represented by FIGURE 2D, and the power of the signal plus noise as 
represented by FIGURE 2A. Then, weights W" can be assigned to each CPPi (i.e., F 
CPPi in FIGURE 2D), based on these estimated signal power and signal plus noise 
power, before forming the integrated signal, ICPP(Ni), and the square root of this 
signal, as depicted in FIGURE 2E. 

Although a preferred embodiment of the method and apparatus of the present 
invention has been illustrated in the accompanying Drawings and described in the 
foregoing Detailed Description, it will be understood that the invention is not limited 
to the embodiment disclosed, but is capable of numerous rearrangements, 
modifications and substitutions without departing from the spirit of the invention as 
set forth and defined by the following claims. In particular, the methods described and 
claimed herein for uplink applications can be used alternatively for downlink 
applications. 
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WH AT IS CLAIMED IS : 

1 . A method for determining, in a radio environment with time dispersion, 
the distance between a radio transmitter and a radio receiver, comprising the steps of: 
transmitting from said radio transmitter M samples of a signal; 
5 receiving in said receiver said M samples together with multipath components 

and noise; 

determining an estimated channel power profile for each of said M samples; 

selecting a first set of N samples from said M samples; 

performing incoherent integration for said estimated channel power profiles for 
10 said first set of N samples to form a first integrated signal; 

if a quality level of said first integrated signal with respect to signal to noise 
is less than a predetermined threshold, selecting another sample from said M samples; 

performing incoherent integration for said estimated channel power profiles for 
said first set of N samples and said another sample to form a second integrated signal; 
15 if a quality level of said second integrated signal with respect to signal to noise 

is greater than or equal to said predetermined threshold, determining a time-of-arrival 
of a maximum level of said second integrated signal; 

entering said time-of-arrival into a time-of-arrival versus frequency of 
occurrence array; 

2 o selecting a second set of N samples from said M samples; 

repeating all of said performing through said entering steps for said second set 

of N samples; and 

determining a minimum value estimated time-of-arrival from said array. 

25 2. The method of Claim 1 , wherein the step of determining a minimum 

value estimated time-of-arrival from said array comprises reading out from said array 
a time value having a first predetermined number of all occurrences being higher time- 
of-arrival values, and a complement of said first predetermined number of all 
occurrences being lower time-of-arrival values. 

3 0 3. The method of Claim 2, wherein said first predetermined number is 

greater than 70%, and said second predetermined number is 100% minus said first 
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predetermined number in percent. 



4. The method of Claim 1 , wherein the receiving step comprises receiving 
said M samples separately in two independent radio receivers to form 2 times said M 
samples. 

5 5. The method of Claim 1 , wherein said array comprises a matrix. 

6. The method of Claim 1, wherein a manual version of said entering step 
comprises entering said time-of- arrival into a time-of- arrival versus frequency of 
occurrence diagram. 

10 7. A method for determining the distance between a mobile station and 

a base station in a mobile network, comprising the steps of: 

transmitting from said mobile station a first plurality of signal samples; 
receiving at said base station said first plurality of signal samples; 
forming an estimated channel power profile value for each of said first plurality 
15 of signal samples; 

selecting a second plurality of samples from each said estimated channel power 
profile value; 

integrating said estimated channel power profile value for said second plurality 
of samples to form a first integrated channel power profile; 
2 0 if a quality level of said first integrated channel power profile is less than a 

predetermined threshold value, selecting at least one additional sample from each said 
estimated channel power profile value; 

integrating said estimated channel power profile for said second plurality of 
samples and said at least one additional sample to form a second integrated channel 
2 5 power profile; 

if a quality level of said second integrated channel power profile is greater than 
or equal to said predetermined threshold value, determining a time-of-arrival for a 
maximum of said second integrated channel power profile; 

entering said time-of-arrival into a time-of-arrival versus frequency of 
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occurrence array; 

selecting a third plurality of samples from each said estimated channel power 
profile value; 

repeating said integrating through said entering steps for said third plurality of 
5 samples; and 

determining a minimum time-of-arrival from said array. 

8. The method of Claim 7, wherein said integrating steps comprise 
incoherent integration. ■ . , 

9. The method of Claim 7, wherein said minimum time-of-arrival 
1 0 represents a shortest distance between said mobile station and said base station. 

10. The method of Claim 7, wherein the step of determining a minimum 
time-of-arrival comprises reading out from said array a time value having a first 
number of all occurrences being higher time-of-arrival values, and a complement of 
said first number of all occurrences being lower time-of-arrival values. 

! 5 11. The method of Claim 1 0, wherein said first number is greater than 0.7, 

and said complement of said first number is 1 minus said first number. 

12. A method for.determining, in a radio environment with time dispersion, 
a distance between a radio transmitter and a radio receiver, comprising the steps of: 
transmitting from said radio transmitter M samples of a signal; 
receiving in said radio receiver said M samples together with multipath 
components and noise; 

selecting a first set of N samples from said M samples to form a mathematical 
function of said N samples, varying the number for N iteratively such that said 
mathematical function of said N samples preserves time information and reduces said 
noise to a first predetermined value, said first predetermined value indicating a desired 
level of said signal relative to said noise; 

creating additional sets of N samples from said M samples, by repeating said 
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selecting step and starting each said repeating step with a different random , and 
independent set of N samples, and terminating said repeating step when a number of 
said additional sets of N samples comprises a second predetermined value; 

determining a time-of-arrival for each of said first set and said additional sets 
5 of N samples; 

storing each said time-of-arrival in a time-of-arrival versus frequency of 
occurrence array; and 

determining an estimated time-of-arrival from said time-of-arrival versus 
frequency of occurrence array. 

10 

13. A method for performing, in a radio environment with time dispersion 
and noise, a determination of an earliest time-of-arrival for a known received signal 
transmitted repetitively as M samples from a remote transmitter, said M samples 
subject to fading such that a time-of-arrival varies between at least two of said M 
15 samples, comprising the steps of: 

' selecting from among said M samples a first set of samples; 
forming a mathematical function from said first set of samples to improve a 
signal-to-noise quality without destroying time-of-arrival information; 

repeating said selecting and forming steps for additional sets of samples in 
2 0 addition to said first set of samples; 

determining a time-of-arrival for each of said first set and said additional sets 
of samples; and 

statistically processing each said time-of-arrival obtained from said 
determining step to obtain an; estimate of said earliest time-of-arrival. 

2 5 14. The method of Claim 13, wherein said first set of samples and said 

additional sets of samples comprise an equal number of N samples. 

15. The method of Claim 13, wherein said first set of samples and said 
additional sets of samples comprise different numbers of samples. 

16. The method of Claim 13, wherein a number of samples used in said 
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first set of samples and said additional sets of samples is selected to obtain a desired 
signal-to-noise quality for each of said first set of samples and said additional sets of 
samples. 

17. The method of Claim 13, wherein a number of samples used in said 
first set of samples and said additional sets of samples is selected to obtain a desired 
signal-to-noise quality defined by h/sigma. 

18. A system for use in determining the distance between a mobile station 
and a base station in a mobile network, comprising: 

a mobile transmitter for transmitting a first plurality of signal samples; 
a base station receiver for receiving said first plurality of signal samples; and 
a processor coupled to said base station receiver, said processor operable to: 
form an estimated channel power profile value for each of said first plurality 
of signal samples; 

select a second plurality of samples from each said estimated channel power 
profile value; t 

integrate said estimated channel power profile value for said second plurality 
of samples to form a first integrated channel power profile; 

select at least one additional sample from each said estimated channel power 
profile value if a quality level of said first integrated channel power profile is less than 
a predetermined threshold value; 

integrate said estimated channel power profile for said second plurality of 
samples and said at least one additional sample to form a second integrated channel 
power profile; 

determine a time-of-arrival for a maximum signal amplitude of said second 
integrated channel power profile if a quality level of said second integrated channel 
power profile is greater than or equal to said predetermined threshold value; 

enter said time-of-arrival into a time-of-arrival versus frequency of occurrence 

array; 

select a third plurality of samples from each said estimated channel power 
profile value; 
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repeat the first integrate operation through the enter operation for said third 
plurality of samples; and 

determine a minimum time-of-arrival from said array. 
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